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We present cluster decay studies of ^^Ni* formed in heavy-ion collisions using different 
Fermi density and nuclear radius parameters proposed by various authors. Our study 
reveals that different technical parameters do not alter the transfer structure of fractional 
yields significantly. The cluster decay half-lives of different clusters lies within ±10% for 
different Fermi density parameters and nuclear radius, therefore, justify the current set 
of parameters used in the literature for the calculations of cluster decay. 
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1 Introduction 



In earlier days, nucleus was considered to have a uniform density and sharp radius. With 
the passage of time, the density distribution was found to be more complicated. Several 
different forms (direct or indirect) exist in literature that can explain these complicated 
nuclear density distributions. The first method is the direct parametrization which in- 
volves the choice of a suitable functional form where parameters are varied to fit the 
experimental data. The two parameter Fermi density distribution is an example of such 
a parametrization. The second method is of indirect parametrization of density distribu- 
tion proceeds via nuclear models. The nuclear models like shell model contains certain 
parameters which are determined by other physical considerations and it is then used to 
calculate the nuclear density distribution without further adjustments. The experimental 
data can be described accurately with two-parameter Fermi density distribution at rela- 
tively low momentum. Among all the density distributions two-parameters Fermi density 
has been quite successful in the low, medium and heavy mass regions. The systematic 
study of charge distributions have been carried out in Refs.[l-3]. We shall use this density 
distribution here. 

Since the nuclear systems obey quantum laws, therefore, their surfaces are not well 
defined. The nuclear density remains constant up to certain distance but fall more rapidly 
close to the surface region where the nucleons are free to move about. The nuclear densities 
provide important information about the structure of nuclear matter at low energies and 
other important information regarding the equation of state at intermediate energies [4, 5]. 

Various methods have been developed for exploring the nuclear structure and radius. 
The electron scattering/ electrically charged particles of high energy are employed as probe 
to explore the proton distribution of the nuclei (i.e charge radii), whereas neutral nuclear 
probes such as neutrons will give the effect of nuclear forces over the nuclear surface 
(i.e. interaction radii). The charge radii are often used to extract the information about 
nuclear radii. The electron scattering experiments shows that the charge distribution 
within a nucleus either follow Fermi trapezoidal shape or modified Gaussian distribution. 
These studies have shown that nuclear charge density does not decrease abruptly but has 
a finite diffuseness. 

A model that uses density distribution such as two parameter Fermi density (as shown 
in Fig. 1) has to rely on the information about nuclear radius (or half density radii Rq), 
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central density po, and surface diffuscncss (a). Interestingly, several different experimen- 
tal as well as theoretical values of these parameters are available in literature [6-11]. In 
addition, several different names such as central radii, equivalent sharp radii, root mean 
square radii etc. have also been used in the literature to define different functional forms. 
The role of different radii was examined in exotic cluster decay half-hves [12] and interest- 
ingly two different forms of radii were found to predict five order of magnitude different 
half-lives within the same theoretical model. Similarly, the use of different values of sur- 
face diffuseness also varies from author to author. The effect of these model ingredients 
on the fusion process at low incident energy have been studied in Ref. [13] and there was 
found that the effect of different radii is more than marginal and therefore this parameter 
should be used with a more fundamental basis. Unfortunately, no systematic study is still 
available in the cluster decay process. In this paper, we plan to study the role of Fermi 
density parameters in the cluster decay of ^^Ni* when formed in heavy-ion collisions. This 
study is still missing in the literature. 

Heavy-ion reactions provide a very good tool to probe the nucleus theoretically. This 
includes low energy fusion process [14], intermediate energy phenomena [15] as well as 
cluster-decay and/or formation of super heavy nuclei [16,17]. In the last one decade, 
several theoretical models have been employed in the literature to estimate the half-life 
times of various exotic cluster decays of radioactive nuclei. These outcome have also been 
compared with experimental data. Among all the models employed preformed chister 
model (PCM) [18-20] is widely used to study the exotic cluster decay. In this model 
the clusters/ fragments are assumed to be pre-born well before the penetration of the 
barrier. This is in contrast to the unified fission models (UFM) [21-23], where only 
barrier penetration probabilities are taken into account. In either of these approach, one 
needs complete knowledge of nuclear radii and densities in the potential. 

Cluster decay of ^^Ni is studied when formed as an excited compound system in heavy- 
ion reactions. Since ^^Ni has negative Q-value (or Qout) and is stable against both fission 
and cluster decay processes. However, if is is produced in heavy-ion reactions depending 
on the incident energy and angular momentum involved, the excited compound system 
could either fission, decay via cluster emissions or results in resonance phenomenon. The 
^^Ni has a negative Qout having different values for various exit channels and hence would 
decay only if it were produced with sufficient compound nucleus excitation energy E^^ {— 
Ecm + Qin)-i to Compensate for negative Qmiti the deformation energy of the fragments 
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Ed, their total kinetic energy (TKE) and the total excitation energy (TXE), in the exit 
channel as: 

EcN =1 Qout I +Ed + TKE + TXE. (1) 

(see Fig. 2, where E^, is neglected because the fragments are considered to be spherical). 
Here Qin adds to the entrance channel kinetic energy Ecm of the incoming nuclei in their 

ground states. 

Section 2 gives some details of the Skyrme energy density model and preformed cluster 
model and its simplification to unified fission model. Our calculations for the decay half- 
life times of ^^Ni compound system and a discussion of the results are presented in Section 
3. Finally, the results are summarized in Section 4. 

2 Model 

2.1 Skyrme Energy Density Model 

In the Skyrme Energy Density Model (SEDM) [7] , the nuclear potential is calculated as a 
difference of energy expectation value E of the colliding nuclei at a finite distance R and 
at complete isolation (i.e. at oo) [7,24]. 

Vn{R)^E{R)-E{oo), (2) 

where E — J H{r)dr, with H(r) as the Skyrme Hamiltonian density which reads as: 

+ ^(^1 + t2)pT + ^(^2 - ti){pnTn + PpTp) 

+^(^2 - 3ii)pVV+ ^isPnPpP 

+ ^(3tl + t2)(PnVVn + PpVVp) 

-^WoipV ■ J + PnV ■ Jn + PpV • Jp). (3) 

Here J = J„ + Jp is the spin density which was generalized by Puri et al. [7], for spin- 
unsaturated nuclei and r = t„ + Tp is the kinetic energy density calculated using Thomas 
Fermi approximation [25,26], which reduces the dependence of energy density H{p,T, J) 

— * 

to be a function of nucleon density p and spin density J only. Here strength of surface 
correction factor is taken to be zero (i.e. A = 0). The remaining term is the nucleon 
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density p = pn + Pp is taken to be well known two-parameter Fermi density. The Coulomb 
effects are neglected in the above energy density functional, but will be added explicitly. 
In Eq. (3), six parameters to, ti, t2, t^, Xq, and Wq are fitted by different authors to obtain 
the best description of the various ground state properties for a large number of nuclei. 
These different parameterizations have been labeled as S, SI, SII, SIII etc. and known 
as Skyrme forces for light and medium colliding nuclei. Other Skyrme forces are able to 
reproduce the data for heavy systems better. The Skyrme force used for the present study 
is SIII with parameters as: to = -1128.75 MeVfm^ ti = 395.00 MeVfm^ t2 = -95.00 
MeVfm^ ts = 14000.00 MeVfm^ xq = 0.45, and Wq = 120.00 MeVfm^ It has been 
shown in previous studies that SIII force reproduces the fusion barrier much better than 
other sets of Skyrme forces for light and medium nuclei. Other Skyrme forces such as 
SKa, SKm, however, are found to be better for heavier masses. 

Prom Eq. (3), one observes that the Hamiltonian density H{p,T,J) can be divided 
into two parts: (i) the spin-independent part Vp{R), and (ii) spin-dependent Vj{R) [7] 
as: 

Vn{R) = J {H{p) - [H^{p^) + H2{p2)]} dr 

j [h{p, J) - Ji) + H2{p2, J2)\ ] dr 

Vp{R) + Vj{R) (4) 



+ 



We apply the standard Fermi mass density distribution for nucleonic density: 

P^^) = ^; R-Ro \ ' - oo < i? < oo (5) 

1 exp I ^^^^ I 

Here po, and "a" are respectively, the average central density, half-density radius and 
the surface diffuseness parameter. The Rq gives the distance where density drops to the 
half of its maximum value and the surface thickness s (= 4.4a) has been defined as the 
distance over which the density drops from 90% to 10% of its maximum value is the 
average central density po- The systematic two parameter Fermi density distribution is 
shown in Fig. 1. 

Another quantity, which is equally important is the r.m.s. radius (r^)^ defined as: 

oo 

(r^>^ = J r'^p{r)dr^in J p{r)r^dh. (6) 
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One can find the half density radius by varying surface diffuseness "a" and keeping r.m.s. 
radius (r^)^ constant or from normalization condition: 

1 



(7) 



The average central density po given by [27] 

3A 



Po = 



1 + 



Using Eq. (5), one can find the density of neutron and proton individually as: 

N Z 



Pn = -JP, 



Pp 



(8) 



(9) 



For the details of the model, reader is referred to Ref. [7]. 

In order to see the effect of different Fermi density parameters on the cluster decay 
half-lives, we choose the following different Fermi density parameters proposed by various 
authors. 

1. H. de Vries et al. [11]: Here, we use the interpolated experimental data [28] of 
Elton and H. de Vries for half density radius i?o and surface thickness a. Using Rq 
and a, central density po can be computed using Eq. (7). This set of parameters is 
labeled as DV. 

2. Ngo-Ngo [6]: In the version of Ngo-Ngo, a simple analytical expression is used for 
nuclear densities instead of Hartree-Fock densities. These densities are taken to be 
of Fermi type and written as: 



Pn,p{R) 



Pn,p(0) 



1 + exp[(i? - C,,p)/0.55] ' 



(10) 



Pn,p(0) are then given by: 



3 iV 1 , , 3 Z 1 

where C represents the central radius of the distribution. 

1 



fill 



C^R 



and 



R 



i?2 

NR^ + ZRp 



(12) 



(13) 
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The sharp radii for proton and neutron are given by, 



with 



= 1.128 fm, ro„ = 1.1375 + 1.875 x 10"^. 
This set of parameters is labeled as Ngo. 



(14) 



(15) 



3. S.A. Moszkwski [8]: The Fermi density parameters due to Moszkwski has central 

density po = 0.16 nucl./fm^, the surface diffuseness parameters a is equal to 0.50 fm 
and radius Rq = 1.15 A^^^. This set of parameters is labeled as SM. 

4. E. Wesolowski [9]: The expressions for Fermi density parameters taken by E. 
Wesolowski reads as: The central density 



Po 



-nRl {1 + {na/Rof} 



(16) 



The surface diffuseness parameters a — 0.39 fm and half density radius, 



Rq — R' 



1 - 



+ 



3\R 



+ 



with 



R' 



1.2- 



0.96 fN-Z 



AV3 \^ A 
This set of parameters is labeled as EW. 



A^/^, and b = -^a. 



;i7) 



(18) 



5. H. Schechter et al. [10]: The value of Fermi density parameters taken by H. 
Schechter et al. can be summarized as: central density po = 0.212/(1 + 2.66A~^/^), 
the surface diffuseness parameters a is equal to 0.54 fm and radius Rq = l.OAA^^^ 
in single folding model for one of the nucleus. This set of parameters is labeled as 
HS. 

In the spirit of proximity force theorem, the spin independent potential Vp{R) of the 
two spherical nuclei, with radii Ci and C2 and whose centers are separated by a distance 
R — s + Ci + C2 is given by 

Vp(R) ^27rR(l)(s), (19) 
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where 

= J {H{p) - [H,{p,) + H2{P2)]} dZ, (20) 

and 

with Siissmann central radius C given in terms of equivalent spherical radius R as 

C = i?-A (22) 

Here the surface diffuseness 6 = 1 fm and nuclear radius R taken as given by various 
authors in the literature [6,29-34]. 

In the original proximity potential [29], the equivalent sharp radii used are 

i? = 1.28^^/^ -0.76 + 0.8/1-^/3 (^23) 

This radius is labeled as Rprox77- 

In the present work, we also used the nuclear radius due to Aage Winther, labeled as 
Raw s-nd read as [30]: 

1.20A^/3 -0.09 fm. (24) 
The newer version of proximity potential uses a different form of nuclear radius [31] 

R = 1.240A^/3 [1 + 1.646^-^ - 0.191A,] fm. (25) 

This radius is labeled as Rproxoo- 

Recently, a newer form of above Eq. (25) with slightly different constants is reported 

[32] 

R = 1.2332A^/=^ + 2.8961^-^/3 - 0.18688^^/^^, fm, (26) 

and is labeled as R^oyer- 

For Ngo and Ngo [6] nuclear radius, we use Eqs. (13)- (15) and is labeled as Rnqo- 
The potential based on the classical analysis of experimental fusion excitation func- 
tions, used the nuclear radius (labeled as Rsass) [33] as: 

R = 1.16/1^/3 - 1.39^-^/3. (27) 

The empirical potential due to Christensen- Winther (CW) uses the same radius form 
(Eq. (27)) having different constants (labeled as Rcw) [34]. 

R = 1.233^^/3 - 0.978/1-^/3. ^28) 
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2.2 The Preformed Cluster Model 

For the cluster decay calculations, we use the Preformed Cluster Model [18-20]. It is 
based on the well known quantum mechanical fragmentation theory [35-38], developed 
for the fission and heavy-ion reactions and used later on for predicting the exotic cluster 
decay [39-41] also. In this theory, we have two dynamical collective coordinates of mass 
and charge asymmetry t] = {Ai — A2)/{Ai + A2) and rjz = {Zi — Z2)/{Zi + Z2). The 
decay half- life T1/2 and decay constant A, in decoupled rj- and i?-motions is 

In 2 

A = = Poi^oP, (29) 
-t 1/2 

where the preformation probability Pq refers to the motion in 77 and the penetrability 
P to i?-motion. The is the assault frequency with which the cluster hits the barrier. 
Thus, in contrast to the unified fission models [21-23], the two fragments in PCM are 
considered to be pre-born at a relative separation co-ordinate R before the penetration 
of the potential barrier with probability Pq. The preformation probability Pq is given by 

Po(A) =1 r ^/bM (^) > = 1 or 2), (30) 

with ip''{r]), 1/ — 0,1, 2, 3, , as the solutions of stationary Schrodinger equation in 77 at 

fixed R, 

riv) = E'-riv), (31) 

solved ai R — Ra — Rmin the minimum configuration i.e. Ra — Rmm (corresponding 
to Vmin) with potential at this i?a- value as V{Ra — Rmin) — Vmin (displayed in Fig. 2). 

The temperature effects are also included here in this model through a Boltzmann-like 
function as 

I m \'= E I ^^(^) I' ' (32) 

i/=0 ^ ^ 

where the nuclear temperature T (in MeV) is related approximately to the excitation 
energy E^^, as: 

E*cj,^Iat'-T, (in MeV). (33) 

The fragmentation potential (or collective potential energy) VR{r]), in Eq. (31) is calcu- 
lated within Strutinsky re-normalization procedure, as 

i=l 

9 



h'' d 



d 



VLDM{Ai, Zi) + 5L/jexp ( -^j 



+ + V,(R), (34) 



where the hquid drop energies {Vldm = B ~ 5U) with B as theoretical binding energy 
of MoUer et ah [42] and the shell correction 5U calculated in the asymmetric two center 
shell model. The additional attraction due to nuclear interaction potential Vjq{R) is 
calculated within SEDM potential using different Fermi density parameters and nuclear 
radii as discussed earlier. The shell corrections are considered to vanish exponentially for 
Eqn ^ 60 MeV, giving T — 1.5 MeV. The mass parameter representing the kinetic 
energy part of the Hamiltonian in Eq. (31) are smooth classical hydrodynamical masses 
of Kroger and Scheid [43] . 

The WKB action integral was solved for the penetrability P [41]. For each 77- value, 
the potential V{R) is calculated by using SEDM for R > Rd, with Rd = Rmin + ^R and 
for R < Rd, it is parameterized simply as a polynomial of degree two in R: 



where Rq is the parent nucleus radius and AR is chosen for smooth matching between 
the real potential and the parameterized potential (with second-order polynomial in R). 
A typical scattering potential, calculated by using Eq. (35) is shown in Fig. 2, with 
tunneling paths and the characteristic quantities also marked. Here, we choose the first 
(inner) turning point Ra at the minimum configuration i.e. Ra = Rmin (corresponding 
to Vmin) with potential at this Ra-value as V{Ra = Rmin) = V min and the outer turning 
point Ri) to give the (5e//-value of the reaction V{Ri,) = Qeff- This means that the 
penetrability P with the de-excitation probability, Wi = cxp{—bEi) taken as unity, can 
be written as P = PiPb, where Pj and P^ are calculated by using WKB approximation, 
as: 



here Ra and Rj, are, respectively, the first and second turning points. This means that 
the tunneling begins at R — Ra {— Rmin) and terminates at R — R^, with V{Rb) — Qeff- 
The integrals of Eqs. (36) and (37) are solved analytically by parameterizing the above 
calculated potential V{R). 




(35) 




(36) 



and 




(37) 
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The assault frequency uq in Eq. (29) is given simply as 




V 



Ro 



(38) 



where E2 — ^Qeff is the kinetic energy of the emitted cluster, with Qeff shared between 
the two fragments and fj, = m( ^'^_^ ) is the reduced mass. 



The PCM can be simplified to UFM, if preformation probability Pq = 1 and the 
penetration path is straight to Qe//- value. 

3 Results and Discussions 

In the following, we see the effect of different Fermi density parameters and nuclear radii 
on the cluster-decay process using the Skyrme energy density formalism within PCM and 



First of all, to see the effect of different Fermi density parameters on the cluster decay 
half-lives, we choose the different Fermi density parameters proposed by various authors 
as discussed earlier. 

Fig. 2 shows the characteristic scattering potential for the cluster decay of ^^Ni* into 
i6q _|_ 40Qg^ channel as an illustrative example. In the exit channel for the compound 
nucleus to decay, the compound nucleus excitation energy E^j^ goes in compensating the 
negative Qout, the total excitation energy TXE and total kinetic energy TKE of the 
two outgoing fragments as the effective Q- value (i.e. TKE = Qeff in the cluster decay 
process). In addition, we plot the penetration paths for PCM and UFM using Skyrme 
force SIII (without surface correction factor, A = 0) with DV Fermi density parameters. 
For PCM, we begin the penetration path at Ra — Rmin with potential at this i?a- value as 
V{Ra — Rmin) — Vmin and ends at it! = Rb, corresponding to V{R — Rj,) — Qeff, whereas 
for UFM, we begin at Ra and end at Rb both corresponding to V{Ra) — V{Rb) — Qeff- 
We have chosen only the case of variable Qe// (as taken in Ref. [44]), for different cluster 
decay products to satisfy the arbitrarily chosen relation Qeff = 0.4(28— | Qout I) MeV, 
as it is more realistic [45]. The scattering potential with SM Fermi density parameters is 
also plotted for comparison. 

Fig. 3(a) and (b) shows the fragmentation potential V{r]) and fractional yield at R = 
Rmin with V{Rmin) — ^ min- The fractional yields are calculated within PCM at T = 
3.0 MeV using various Fermi density parameters for ^^Ni*. From figure, we observe that 



UFM. 
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different parameters have minimal role in the fractional mass distribution yield. The fine 
structure is not at all disturbed for different sets of Fermi density parameters. 

We have also calculated the half-life times (or decay constants) of ^^Ni* within PCM 
and UFM for clusters >^^0. For ^^O, the cluster decay constant varies by an order of 
magnitude ten. The variation is much more with SM parameters. In the case of UFM, 
variation is almost constant. 

In Fig. 4, we display the cluster decay half-lives logTi/2 for various Fermi density 
parameters using PCM. There is smooth variation in half-life times with all the density 
parameters except for SM parameter. The trends in the variation of cluster half-life times 
(or decay constants) are similar in both PCM and UFM, but in case of UFM decay 
constants are more by an order of ten. In SM the decay constants are larger by an order 
of 14. 

In order to quantify the results, we have also calculated the percentage variation in 
logTi/a as: 

[log Ty,] % = ilogT,/,)^y ^^^^ 

where i stands for the half-life times calculated using different Fermi density parameters. 
The variation in the cluster decay half-lives is studied with respect to DV parameters. In 
Fig. 5(a) and (b), we display the percentage variation in the half-life times within both 
the PCM and UFM models as a function of cluster mass A2 using Eq. (39). For the 
PCM these variation hes within ±5% excluding SM parameters, whereas including SM 
parameters it lies within ±13%. In the case of UFM half-lives hes within ±1.5% for all 
density parameters except of SM. For SM parameters variations lie within ±9%. 

Finally, it would be of interest to see how different forms of nuclear radii as discussed 
earlier would affect the cluter decay half-lives. 

In Fig. 6, we display the characteristic scattering potential for the chistcr decay of ^^Ni* 
into ^^Si + ^^Si channel for Rsass and Ruoyer forms of nuclear radius. In the exit channel 
for the compound nucleus to decay, the compound nucleus excitation energy Eqj^ goes in 
compensating the negative Qout, the total excitation energy TXE and total kinetic energy 
TKE of the two outgoing fragments as the effective Q-value. We plot the penetration 
path for PCM using Skyrme force SIII (without surface correction factor, A = 0) with 
nuclear radius R^ass- Here again, we begin the penetration path at Ra = Rmin with 
potential at this i?a-value as V{Ra = Rmin) = ^mm and ends at R = Rb, corresponding 
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to V{R = Rb) = Qeff for PCM. The Qeff are same as discussed earlier. 

Fig. 7(a) and (b) show the fragmentation potentials V{r]) and fractional yields at 
R — Rmin with V{Rmin) — V min- The fractional yields are calculated within PCM at 
T — 3.0 MeV for ^^Ni* using various forms of nuclear radii. Prom figure, we observe 
that different radii gives approximately similar behavior, however small changes in the 
fractional mass distribution yields are observed. The fine structure is not at all disturbed 
for different radius values. 

We have also calculated the half-hfe times (or decay constants) of ^^Ni* within PCM 
for clusters >^^0. The cluster decay constant for nuclear radius due to Bass varies by an 
order of magnitude 10^, where as order of magnitude is same for other radii. In Fig. 8, 
we display the cluster decay half- lives logTi/2 for various nuclear radii taken by different 
authors as explained earlier using PCM. One can observe small variations in half-life 
times. 

In order to quantify the results, we have also calculated the percentage variation in 
log 71/2 as: 

[iogr,j./o^ "°«^-;;>--"°e^''''""" xi°°. m 

L J (logTi/2)^'«°s'«'' 

where i stands for the half-hfe times calculated using different forms of nuclear radii. The 
variation in the cluster decay half-lives is studied with respect to radius formula given 
by Royer Rijoyer- In Fig- 9, we display the percentage variation in the half-life times for 
PCM as a function of cluster mass A2 using Eq. (40). These variation lies within ±7% 
excluding Bass radius where it lies within ±10%. 

4 Summary 

We here reported the role of various model ingredients as well as radii in the cluster decay 
constant calculations. Our studies revealed that the effect of different density and nuclear 
radius parameters on the cluster decay half-life times is about 10%. Our study justify the 
use of current set of parameters for radius as the effect of different prescriptions is very 
small. 
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Figure 1: The systematics diagram for two parameter Fermi density. 
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Figure 3: (a) The fragmentation potential V{r)) and (b) calculated fission mass distribu- 
tion yield with different density parameters at T = 3.0 MeV. 
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Figure 4: The variation of logTi/2 (sec) using different density parameters for PCM. 
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Figure 7: Same as fig. 3, but for different radii. 
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Figure 8: Same as fig. 4, but for different radii. 
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Figure 9: Percentage variation of logTi/2 for different forms of radii with PCM only. 
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